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I. INTRODUCTION 



Resonances are one of most interesting phenomena in scattering of quantum particles. 
0> ' The problem of definition and studying resonances is payed a lot of attention both in physical 
and mathematical literature (see, e. g., the books [0 ^ [§[)• Main difficulties connected with 
a rigorous definition of resonance are explicitly emphasized by B. Simon in his survey 0. 
O I A presence of such difficulties is obliged, first of all, to the fact that in a contrast to the 
usual spectrum, resonances are not an unitary invariant of an operator (Hamiltonian of a 
quantum system). The generally accepted interpretation of resonance as a complex pole 
- of the scattering matrix continued analytically on unphysical sheet(s) of the energy plane, 
', goes back to the known paper by G. Gamow For radially symmetric potentials, such 
I an interpretation of the two-body resonances has been rigorously approved by R. Jost [IT . 
Beginning from E.G. Titchmarsh |]12| resonances are considered as well as poles of analytic 
continuation of the Green function (or its matrix elements between suitable states 0, |^). 
A survey of different physical approaches to studying three-body resonances may be found 
e. g., in 1^ and |T^. 



At the moment, one of the most effective approaches to practical calculation of resonances 
is the complex scaling method |T3[ (see also 0). This method is applicable to the 
few-body problem in the case where interaction potentials between particles are analytic 
functions of coordinates. The complex scaling gives a possibility to rotate the continuous 
spectrum of the A^-body Hamiltonian in such a way that certain sectors become accessible for 
observation on unphysical sheets neighboring with the physical one. At the same time, the 
real discrete spectrum of the Hamiltonian stays fixed during all the scaling transformation. 
Resonances in the sectors above turn out to be extra discrete eigenvalues of the scaled 
Hamiltonian Thereby, when searching for resonances one may use standard methods to 
find discrete spectrum. Practical applications of the complex scaling method to concrete 
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problems may be found, in particular, in the recent papers [0 - |T^. Alongside with the 
complex scaling, another methods are used for calculations of three-body resonances which 
are based in particular on solving the momentum space Faddeev integral equations 



19| continued through the cut (see the survey [jl3| by K. MoUer and Yu.V. Orlov and 



the literature cited therein). In this approach, resonances are searched for as poles of the 
T-matrix. 

The present paper is devoted to developing a method to calculate three-body resonances 
using the recently found explicit representations disclosing a structure of the T-matrix on 
unphysical sheets as well as analogous representations for the scattering matrix and re- 



solvent ||2y], These representations were obtained in supposition that the interaction 
potentials were pairwise and falling-off in the coordinate representation not slower than 
exponentially. According to the representations ||2T[, the matrix Nl{z) = {Mai3{z)}, 



a,P = 1,2,3, constructed of the operator T{z) Faddeev components |T8|, |T9l, is explicitly 
expressed on unphysical sheet 11; of the energy z plane in terms of this matrix itself taken on 
the physical sheet and a certain truncation Si{z) of the total three-body scattering matrix 
S{z). Character of the truncation is determined by the index (number) / of the unphysical 
sheet concerned. Respective representations for analytic continuation of the matrix S{z) 
and resolvent R{z) follow immediately from the representations for M(z)|jj^. A main con- 
sequence of the representations admitting direct practical applications, is the fact that the 
T-matrix and resolvent as well as the scattering matrix have nontrivial singularities on un- 
physical sheet 11; exactly at those values of the energy z where the corresponding matrix 
Si{z) has zero as eigenvalue. It is important that Si{z) is considered on the physical sheet 
only. Therefore, one can provide a search for resonances (poles of M(2;), S{z) and R{z) ) 
on a certain unphysical sheet 11; keeping z always on the physical one and calculating only 
a position of zeros of the operator- valued function Si{z). For all this, one can use any 
method allowing to calculate (on the physical sheet) amplitudes of the processes necessary 
to construct the truncation Si{z). 

In the present paper, the matrices Si{z) are computed on the base of the numerical 
algorithm [^] elaborated to solve the Faddeev differential equations in configuration space 
(see the book |jl9|], survey and references therein). Certainly, when computing the 
amplitudes on the physical sheet one has to extend the Faddeev differential formulation of 
the scattering problem as well on the complex values of z. It should be noted that, in the 
holomorphy domain (see [^) of the amplitudes, the differential formulation stays to be 
correct. 



Unfortunately, the algorithm |T9|, [^, p3[ (see also p4| - [^) has been worked out in 



details only for the processes (2 — > 2,3). Thus, there may be computed in practice only 
the amplitudes of elastic scattering and rearrangement for the processes (2 — > 2) and the 
breakup amplitude into three particles. A knowledge of these amplitudes is sufficient to 
compute those truncations Si{z) of the three-body scattering matrix S{z), zeros of which 
are "responsible" for resonances situated on the so-called two-body unphysical sheets, i. e. 
those sheets of the energy Riemann surface where the parameter z may be guided going 
around the pair thresholds only. As a concrete application of the method concerned we 
make a search for resonances in the nnp system and in a model system of three bosons with 
the nucleon masses. 

Let us describe shortly structure of the paper. 
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In Sec. H, we introduce main notations and formulate the explicit representations ||2T| for 



the unphysical-sheet three-body T-matrix, scattering matrix and resolvent which are used 
then to approve the numerical method of the work. 

In Sec. m the system nnp and a three-boson system are considered. Formulations of 
the boundary- value problems [|19|, ||2^, are given for the Faddeev partial differential 



equations, corresponding to the processes (2 — > 2, 3) in these systems and going out to a 
domain of complex energy values in the physical sheet. Numerical method to solve these 
problems is described. Truncated (partial) scattering matrices are constructed in terms of 
the amplitudes for elastic scattering (2 — > 2). Zeros of these matrices represent resonances 
(including virtual poles) on the unphysical sheet connected with the physical one by crossing 
the continuous spectrum interval between the deuteron energy and three-body threshold. 
Results of numerical computations are exposed. 

We conclude this introduction with some notation. Throughout the paper, we understand 
by y/z - A, z eC, XeR, the main branch of the function {z — A)^/^. By p we denote the 
unit vector in the direction p G M", p = ^, and by S"'\ the unit sphere in M^ p G S-\ 

II. EXPLICIT REPRESENTATIONS FOR -MATRIX, SCATTERING MATRIX 
AND RESOLVENT ON UNPHYSICAL SHEETS 

The scattering matrix, -matrix and Hamiltonian resolvent (Green function) for a 
quantum-mechanical system are stringently connected with each other. Therefore, all these 
three objects, considered as functions of energy, have usually the same Riemann surface. 
Such a fact takes place at least in the multichannel scattering problem with binary channels 



and in the three-body problem with quickly decreasing interactions I^H], |2T[]. In a "first 
approximation", the structure of the Riemann surfaces in these problems coincides. The 
thing is that the branching points (in real axis) are stipulated to these surfaces in the both 
problems by the same reason, namely by a presence of the Cauchy-type integrals in the 
Lippmann-Schwinger or Faddeev equations. In the equations considered in the momentum 

representation, the Cauchy-type integrals are engendered by the kernels ^ ^ ^ with A, 

the channel thresholds and p, p G M", the respective channel momentum variables. In the 
case of the channels (2 — > 2, 3) in three-body problem and the odd- dimensional channels 
(i. e., with n odd) in the matrix multichannel problem, the thresholds A turn out to be the 
branching points of the second order. Even- dimensional channels in the multichannel prob- 
lem as well as the channel (3 — >■ 2, 3) in the three-body problem give logarithmic branching 
points (see 0). 



The method used for calculation of resonances in the present work, is based on the 
explicit representations |^ for analytic continuation of the T-matrix, scattering matrices 



and Green function on unphysical sheets keeping true at least for a part of the three- 
body Riemann surface. To describe this part we introduce the auxiliary vector-function 

/(^) = {fo{z))Jl,l{z),...Jl,ni{z),f2,l{z),...j2,n2{z),f3,l{^)^-j3,n3{z)) with fo{z) = \n Z 

and fa,j{z) = {z—XajY^"^- Here, by Xaj we understand respective bound-state energies of pair 
subsystems a, a = 1, 2, 3, j = 1, 2, Ua, numerated taking into account their multiphcities. 
It is supposed that < oo. The sheets 11/ of the vector-function f{z) Riemann surface 3fJ 
are numerated via multi-index / = (Zq, /i.m, ^2,1, ^2,n2! ^3,i; ^a.ns), where = 
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if the sheet 11; corresponds to the main (arithmetic) branch of the square root {z — A^)^/^. 
Otherwise, laj = 1 is assumed. Value of Iq coincides with the number of branch of the 
function In In z = In + i 27r/o + ^0 where = arg z. For the physical sheet identified by 
^0 = la,j = 0, a = 1, 2, 3, j = 1, 2, Ua, we use the notation IIo. A "sticking" of the sheets 
ri; of the surface 3? is realized in intervals between neighboring thresholds along rims of the 
cut along the continuous spectrum. A detailed description of the surface 3fJ may be found 



in 28 



The type 3? surfaces without extra branching points arise in the multichannel problems 
with binary channels only. Structure of the total three-body Riemann surface is essentially 
more complicated. For instance, the sheets 11/ with Iq = ±1 have additional branching points 
corresponding to resonances of the two-body subsystems. In the sheets 11/ with Iq = 0, one 
can discover (in the left half-plane) logarithmic branching points of a kinematical origin. The 
part 3?^^-* of the total three-body Riemann surface where the representations are valid, 
consists of the surface 3? sheets 11/ identified hj Iq = (such unphysical sheets are called 
two-body sheets) and two three-body sheets identified by /q = ±1 and laj = I, a = 1,2,3, 
j = 1, 2, ...,na- Note that 3?'-^'' includes all the unphysical sheets neighboring with the physical 
sheet Ho. 

In the case of the nnp-system and a model three-boson system considered below, the 
surface SfJ^^-* is shown schematically in Fig. |l|. In this case a single pair threshold, the 
deuteron energy A = i^^ is present only. Therefore, the index / of the sheets 11/ consists of 
two components only: / = (IqJi). In the terminology accepted, the sheet 11(0,1) (see Fig. |I]) 
is a two-body sheet, but the sheets II(_i i) and 11(1 1) are three-body ones. 

Construction of the representations for the T-matrix consists of the following stages. At 
the first step, one carries out analytic continuation on unphysical sheets, of the absolute 
terms and kernels of the Faddeev integral equations for the components Maf^i^z) (the con- 
tinuation is understood in the sense of distributions). As the absolute terms as the kernels 
after continuation are expressed in terms of the pair T-matrices and scattering matrices 
taken on the physical sheet. Transforming the Faddeev equations continued, one finds the 
kernels Mai3{P, P' , z)\^^jj^ can be explicitly expressed in terms of these kernels themselves 
taken on the physical sheet IIq in their off-shell and/or half-on-shell variants. It is supposed 
in the last case that the first argument P of the kernels Mai3{P, P', taken on the 

energy ("mass") shells |Pp = z or \pa\'^ = z — X^j, j = 1,2, ...,na- Here, we use the nota- 



tions P = {ka,Pa} with ka,Pa, tt = 1, 2, 3, Standing for the standard relative momenta |]T9 



Transferring in the expressions obtained all the off-shell terms to the l.h. part and inverting 
an operator arising there, one comes to a closed system of equations for the half-on-shell 
components Mq,^(P, P', z)\^^jj^. This system admits an explicit solution using the terms of 
the physical sheet only. As a result one gets the following representation^ for the matrix 
Nl{z) = {Ma/siz)} ,a, /3 = 1,2, 3, continued on the sheet 11/: 



For the sake of simplicity, we write here these representations as well as representations for the 
scattering matrix and resolvent (see below Eqs. (^) and (^), respectively) for the case of spinless 
particles only. A direct generalization of the representations on the case of spin particles causes no 
difficulties. 
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FIG. 1. Physical, Ho, and neighboring unphysical, H^g,!), x) and n(-_i i), sheets of the 
Riemann surface in the three-nucleon (nnp) and model three-boson problems. The near (in figure) 
rim of the cut in the sheet Ho is identified in the interval between the deuteron energy and 
three-body threshold z = with the remote rim of the cut in the sheet 11(0,1) ■ On the contrary, 
the remote rim in the sheet Ho is identified in this interval with the near rim in the sheet H^q,!)- 
In the interval (0, +oo), one identifies the near rim of the cut in the sheet Hq with the remote rim 
in the sheet n^j^ j^-j. Respectively, the remote rim in the sheet XIq is identified in this interval with 
the near rim in the sheet n(_i i). On the sheet 11(0,1); the virtual level Zres(^H) in the s-state of 
the nnp-system is marked, and the resonance Zres(3B) in the model three-boson system. Also, the 
complex conjugate pole Zres(3B) is marked. 
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M(z)|n, = M(;.) - B^{z)A{z)LSr\z)LB{z). (1) 
Here, the factor A{z) is the diagonal matrix, 

A{z) = dmg{Ao{z),Ai^i{z),...,Ai^n.^{z)}, 



combined of the functions Ao{z) = —niz^ and A^j = —T^i^z — X^j. Notations L and L are 
used for diagonal number matrices whose nontrivial elements are the sheet 11/ indices: 

L = diag{/o, h,!, /s.ng} and L = diag{|/o|, hns}- 
By Si{z) we understand a truncation of the three-body scattering matrix: S{z), S{z) : Q —> 

" ^ 3 nc 

g, g = L2{S^) © © ^2(5*^), defined by the equation 

a=l j=l 

Si{z)=i + L[S{z)-i]L 
where I is the identity operator in g. Also, we use the notations 

B(^) = ( + v] ) = {m{z)Q^jI [v + MT]M/Jt) . 

Here, v = diag{wi, V2, f 3} with Va, the pair potentials, a = 1, 2, 3. At the same time, 

1 1 

1] = (1, 1, 1), T = I 1 1 I and ^ = diag{^i,*2,^3} 

1 1 

where "^a, a = 1, 2, 3, are operators acting on / = (/i, /2, /n ) G © 1^2(1^^) as 

i=i 

where, in turn, ipaj is the bound state wave function of the pair subsystem a corresponding 
to the level Xaj- By we denote operator adjoint to Notation Jo{z) is used for operator 
restricting a function on the energy-shell |Pp = z. The diagonal matrix- valued function 

3i{z) = diag{Ji,i(2;), J3,„3(2;)}, 

consists of the operators Ja,j{z) of restriction on the energy surfaces = z — Xaj- The 
operators fl\ ^o{z) and j\{z) represent the "transposed" matrices Q, Jo(-z) and Ji{z), re- 
spectively. Operators Jo(^) and j\{z) act in the expression for (as if) to the left. 

Representations for the scattering matrix and resolvent on unphysical sheets are an 
immediate consequence of the representations ([l|) for the matrix M(2;)|]-[^. 

With details omitted (see |PT|) the representations for the scattering matrix S{z) read 

S{z)\^^ = £{l) {/ + Sr\z)[S{z) - i]e{l)]£{l). (2) 

Here, £ = diag{£o5 £1,1, •••5 ^^s.na} where £0 is the identity operator in L2{S^) if /q = and £q, 
the inversion, (£^o/)(-P) = f{—P), if ^0 = ±1- Analogously, £a,j is the identity operator in 
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L2{S'^) for laj = and inversion for laj = 1. Notation e{l) is used for diagonal number matrix 
e(/) = diagjeo, ei^i, es^ng} with nontrivial elements Caj = 1 if laj = and Caj = —1 if 
la J = 1; for all the cases Cq = 1. 

Analytic continuation R{z)\jj^ of the resolvent R{z) = {H — z)^^ of the three-body 
Hamiltonian H admits the representation 

R(z)\^^ = R + Q^ALSf'LQ. (3) 

Here, 

« = ( J.*.|f- TMlpt ) . Q* = (1^ - . "1^ - ft=MT|*jt) 

with 1^ = ^1-1-^2 + ^3 and -Ro(^) = {Hq — z)~^, the resolvent of the kinetic energy operator 
Hq for the system under consideration. 

There were holomorphy domains n['^°''' found in for the truncated scattering matrices 
Si{z) in the physical sheet. Representations (|I|) and are valid in the same domains. 

It follows from the representations (|l]) — that the resonances (the nontrivial poles of 
M(2;)|j-[^ , S{z)\jj^ and R{z)\jj^) situated on the unphysical sheet 11; are in fact those points 
z = z,:es on the physical sheet where the matrix Si{z) has zero as eigenvalue. Therefore, a 
calculation of resonances on the unphysical sheet Hi is reduced to a search for zeros of the 
truncation Si{z) of the total three-body scattering matrix S{z) in the physical sheet. At the 
numerical search for the resonances, one can use any method allowing to find analytic contin- 
uation on the physical sheet, of the elastic scattering, rearrangement or breakup amplitudes 
necessary to construct the respective truncation Si{z). 



III. RESULTS OF SEARCH FOR RESONANCES IN THREE-NUCLEON (nnp) 

AND THREE-BOSON SYSTEMS 

One of the most effective methods for a study of concrete three-particle systems is the 



numerical algorithm [^] (see also [Q, |]T9[ and references therein) based on the Faddeev 
differential equations for the wave function components in the configuration space. This 
method gives opportunity of comparatively easy calculations of the scattering wave func- 
tions and respective amplitudes for (2 2, 3) processes. An extension of the differential 
formulation to a domain of complex energies enables us to calculate analytic continuation of 
the (2 — > 2, 3) amplitudes on the physical sheet. This is quite enough for us to construct 
the truncated s-state scattering matrices Si{z) whose zeros are resonances on the respective 
two-body unphysical sheets 11/ with /q = 0. 

As to a search for resonances on the three-body unphysical sheets 11; with Iq ^ (in 
particular on the sheets n(_i i) and 11(1,1) in the three-nucleon problem), the situation is 
much more complicated. To construct Si{z) in this case, one has, alongside with (2 — > 
2), to calculate amplitudes for the processes (3 — > 2,3) with three asymptotically free 
particles in an initial state. Unfortunately, the reliable practical methods for calculation of 
the processes (3 — > 2, 3) are not developed so far even for the real energies. In particular, 
when using the differential formulation of the scattering problem, one has from the very 
beginning to separate explicitly contributions to the Faddeev components not only from 
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single-rescattering but also from double-rescattering processes [O]. In the last case, one 



has to take into account explicitly (see [19|) a presence of the "light" and "shadow" zones 
for the correspondent waves, and to use the Fresnel integral for description of intermediate 
regimes. Moreover, computations of the (3 — > 2, 3) amplitudes have to be carried out for 
many different directions of the incident momentum P. This circumstance is itself a factor 
enlarging numerical complexity as compared with the case of the processes (2 — > 2, 3). 

So, in the present work we restrict ourselves with a search for the nnp system resonances 
situated on the unphysical sheet 11(0,1) only, connected with the physical one by crossing 
the continuum spectrum interval [E^, 0) between the deuteron energy z = E(i and breakup 
threshold z = (see Fig. |1|). To construct the truncated scattering matrix 5'(o,i)(-2) (in 
accordance with ([l|)-(|^) just its zeros represent resonances on the sheet 11(0,1)), it suffices to 
calculate the elastic nd scattering amplitudes only. 

To find these amplitudes we use the two-dimensional Faddeev integro-differential equa- 
tions \W\ being a result of the partial and angular analysis of the Faddeev differential 



equations. In addition we make an assumption, rather usual in nuclear physics, that the 
nucleons involved interact in the s-state only. With this assumption the partial equations 
become exact. The Faddeev component Ul corresponding to the total spin S = 3/2 and 
total orbital momentum L, satisfies the equation 

{HL-z)Ul{x,y) = Vt{x)^l{x,y) (4) 
q2 q2 j^n _|_ i\ 

with Hl = — TT-pT — TTo "I 9 5 the partial Laplacian and x,y, the absolute values of 

the Jacobi vectors [|l^. Respective partial component \l/^(x, y) of the total wave function 
for the nnp-system is expressed by the function Ul{x,y) as 

1 

<ill{x, y) = Ulix, y)-^Jdu h\x, y, u) Ul{x', y') (5) 

-1 



where x' = y ^ + ^ 1/^ — xyu and y' = y ^ + ^ 2/^ + xyu. The geometric func- 
tion is given by 



= Etttt^^^W fv^cos^)'(sinr^-^ 

^x'y'\2sme') k\{L - k)\ ^ ^\ J ^ 

y y' 

with Pkiu), the Legendre polynomial of the order k and 9 = arctg— , 9' = arctg — . The 

X x' 

factor Vt[x) represents a triplet part of the nucleon-nucleon interaction potential. 
Function U'l{x,y) satisfies the boundary conditions 

f^l(^,2/)Uo = and f/l(x,y) 1^^0 = 0- (6) 
It satisfies also the asymptotical conditions 
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p >oo 



+Al{z 



3l{,\/z- EdV) + al{z) exp <; i 

exp \ i^p + 
,0) ^ ^ 



\^ z - EdV + i 



ttL 



+ 



(7) 



where ji is the Bessel spherical function of the order L and p = a/x^ + |/^. By ipdix) 
we denote the deuteron wave function and by a'Kz), the quartet amphtude for elastic nd- 
scattering in a state with angular momentum L. The function A\{z, 9) represents the partial 
Faddeev component of the breakup amplitude of this system into three particles. 

The Faddeev integro-differential equations for the doublet (S = 1/2) nrf-scattering read 

as 



{Hl - z)Ut(2),L{^,y) = Vti^s){x)'^U2),L{^,y) 



(8) 



where , \E'^ 



is expressed by the vector 



as 



uf + 



duh^{u)BUl, B 



( 1 

3 



v 



3 

"3 



1 
5/ 



(9) 



By Vs we understand a singlet part of the nucleon-nucleon interaction potential. 
The Faddeev partial components Ufj^ satisfy as [/|, the boundary conditions 



1=0 



0, utA^.y) 



J/=0 



0, 



(10) 



and have the asymptotics 



p ^oo 



jdV^- Edy) + aiiz) exp \ i^/z - E^y + 



exp ^iy/zp + 



VP 



1,2. 



(11) 



Here, a'l^z) is the doublet elastic nrf-scattering amplitude in the state of the nnp-system 
with angular momentum L. At the same time, the functions Af and A21 represent the 
partial Faddeev components of the breakup amplitude for this state. Remind that the 
physical breakup amplitudes in quartet and doublet states are expressed by the amplitudes 
respectively, A^ and Af j^, A^^ via relations analogous to and (|^) (see e. g., |T9|). 

A component of the truncated scattering matrix >S'(o,i) in the state with a fixed angular 
momentum L is diagonal. Its nontrivial quartet, s\{z), and doublet, sf^(z), elements are 
given by 



sl{z) 



1 + 2za^(z). 



\{z) = l + 2ia'i{z). 
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TABLE I. Convergence of the nnp-system virtual level z^esi^H) 



N 






^(3) 


Hit -L-Li^ 


n f TTI 
Hmax 1 -LJ-Ai 


z (^H) MpV) 


40 


10 


20 


30 


0.40 


19.0 


2.9296 


60 


15 


30 


45 


0.30 


23.9 


2.8229 


80 


20 


40 


60 


0.25 


29.9 


2.7565 


120 


30 


60 


90 


0.15 


34.0 


2.7434 


160 


40 


80 


120 


0.10 


39.1 


2.7282 


180 


45 


90 


135 


0.08 


39.5 


2.7275 



Therefore, the sheet 11(01) resonances in the state L of the n?7,p-system are in fact, zeros of 



the sca/ar functions and s'Kz). 

When solving the boundary- value problems (H), (||), (|^) and (H), ( PI) numerically 
at complex energies z, we use the same algorithm ||2^, as at real z. First, we 
make a finite-difference approximation of the problems above in polar coordinates p, 6. As 
the grid points we take the intersection points of the arcs p = pi, i = 1,2,..., Np, and rays 
9 = 9j, j = l,2,...,Ng, 9j < 9j+i. At po = and given pi (see Table |), the grid points 
in p for 2 > 2 are chosen in such a way that p^+i = pj + a(pj — Pi_i) with a parameter a 
(acceleration) not depending on the number of a grid point. All the results exposed below, 
are related to a = 1,01. The interval [0,7r/2] where the variable 9 changes, is divided by 

(k) 

special points with numbers iV^ , = 1, 2, 3, into three subintervals, inside which the grid 
points in 9 are distributed uniformly. When going from one interval to another (in direction 
of the parameter 9 rise), the grid step is divided by two. The choice of the grid described, is 
explained by a necessity to take into account essentially more quick change of the Faddeev 
component values in the domain where p and/or x [1^ are small. Usually, we chose the 
numbers of points in 9 and p the same, Ng = Np. A maximal value of the parameters Nq, 
Np has equaled to 180. With these values, the cut-off radius pmax = PNp has reached 39 fm. 
Typical values of the grid parameters are given in Table | where N = Ng = Np. 

As a nucleon-nucleon interaction, the Malfiiet-Tjon potential MT I-III is used in its 



initial version [29 



Having solved the problems (0), (|^), (|^) and (||), ([lOD, (pl|) we calculate the functions 
s\{z) and s'l^z). Resonances, considered as roots of these functions in the complex plane, 
are found using the Newton method with a three-point approximation of derivative. 

As a test of the computer code we have calculated the bound-state energy of the 
nucleus as a pole of the function s'Kz). More precisely, this pole was found as a root of 
the inverse amplitude l/a'l^z). Beginning from the grid dimension 80x80, we obtained for 
Et the value —8.55 MeV. Hereafter all the energies are given with respect to the breakup 
threshold. Note that the value stated for Et is in a good agreement with known results on 
Et in the MT I-III model (see pj). 

When searching for the nnp resonances on the sheet H(o,i) at L = 0, 1, we have inspected 
a domain of a range about 10 MeV in vicinity of the segment [E^, 0] in the complex z plane. 
Especially carefully we studied a vicinity of the points z = —1.5 ± 0.3 -|- i(0.6 ± 0.3) MeV 
interpreted in the recent works ||3^, |3ll] as a location of an exited state energy of the 
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FIG. 2. Surface of the function |so(-z)|, so{z) = Sq^{z), in the model system of three bosons 
with the nucleon masses. The potential V'-^{r) is used with the barrier VJ, = 1.5 MeV. Position of 
the resonance 2^res(3B) corresponds to the minimal (zero) value of |5q"^(z)|. 

nucleus. Unfortunately we have succeeded to find only one root z^esi^H) of the function 
Sq{z), corresponding to the known virtual state of the nnp system at the total spin S=l/2. 
Position of this root for different grids is shown in Table H illustrating a degree of convergence 
in the method used. As one can see from the table, for the maximal of the grids examined, 
the grid 180 x 180, we have found z,-es{^H) = —2.728 MeV. This means that the calculated 
virtual level is situated 0.504 MeV to the left from the nd threshold Ed = -2.224 MeV (in 
the MT I-III model p9[). Note that the shift Ed — z^csi^^H) found from experimental data 



on nd scattering, is 0.515 MeV (see [T^). Its value computed in a separabilized MT I-III 
model on the base of the momentum space Faddeev equations, equals to 0.502 MeV [|1^. As 
to the resonance 2 = — 1.5 ±0.3 + ^(0. 6 ±0.3) MeV at L = discussed in it is quite 

possible to be situated not on the two-body sheet 11(0,1) but on the three-body sheets n(_i,i) 
or 11(1,1) (see Fig. |l]). Due to the reasons mentioned at the beginning of the section, we may 
unfortunately neither approve nor disprove this hypothesis. As it should be expected (see 
the data on three- nucleon resonances in 0]), we have failed to find any resonances in the 
quartet state at L = as well as at L = 1. 

Also, we have studied a model three-body system including identical spin-zero bosons 
with masses of the nucleon. Remind that the Faddeev integro-differential equations for 
such a three-boson system read exactly as the equations (^ except a necessity to replace 
the factor —1/2 in the expression (|^) with unity. The boundary conditions for respective 
partial Faddeev components Uf^{x,y) have the form (D and (0) where instead of a^z) and 
A'^j^{9,z), one has to substitute a^i^{z) and A^{9,z). 
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Im2 , MeV 




5.5 



Re z, MeV 

FIG. 3. Trajectory of the resonance 2;rcs(3B) on the sheet 11(0,1) iii the model system of three 
bosons with the nucleon masses. The potential V^{r^ is used. Values of the barrier VJ, in MeV are 
given near the points marked on the curve. 
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Component of the truncated scattering matrix 5'(o,i) for the three-boson system is given 
in the state with the angular momentum L by 

sf{z) = l + 2^af{z). 

Resonances on the sheet 11(0,1) in this state are roots of the equation sf^i^z) = considered 
on the physical sheet. 

In the three-boson problem we restrict ourselves with a treating the s-state only and 
thereby, with a searching for zeros of the function s^{z) at L = 0. As a pairwise interaction 
between the bosons we have used the Gauss-type potential supplied with an additional Gauss 
repulsive barrier term, 

V'^{r) = Voexp[-/ior2] + VJ,exp[-/ib(r - rtf] 

where the values Vq = —55 MeV, /xq = 0.2 fm~^, fih = 0.01 fm~^, r;, = 5 fm have been 
fixed and the barrier amplitude Vb varied. A resonance (with non-zero imaginary part) on 
the sheet 11(0,1) arises in the system concerned just due to the presence of the barrier term. 
Example of a surface of the Sq^{z) absolute value for the barrier amplitude Vb = 1.5 MeV 
is shown in Fig. ^ (for a 80x80 grid). A trajectory of the resonance 2;res(3B) (a zero of 
the function sp(2)) is shown for the changing barrier Vb in Fig. ^ This trajectory was 
watched for the barrier Vb decreasing in the interval between 1.5 MeV and 0.85 MeV. When 
drawing the trajectory, we have used a 160 x 160 grid. It can be seen from Fig. ^ that the 
behavior of the resonance Zres(3B) is rather expected: with monotonously decreasing real 
part, the imaginary part of the resonance changes also monotonously. In Fig. ^ we plot both 
the trajectories of the resonance real part Re2res(3B) and two-boson binding energy E^. As 
one can see from this figure, the value of | Re2;res(3B)| increases more quickly than \Ed\, 
coinciding with \Ed\ at H = 0.97 MeV. 

Trajectory of the resonance concerned in the lower complex half-plane is symmetric to 
the curve shown in Fig. ^ with respect to the real axis. Respective points, symmetric to 
those marked in Fig. |^, correspond to the same values of Vb- For Vb ~ 0.95 MeV the energy 
2;res(3B) becomes real. In this point the resonances 2;rcs(3B) and Zres(3B) "collide" turning 
for Vb < 0.95 MeV into a couple of virtual levels. 
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